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ABSTRACT 



An experiment was conducted in the Florida Keys from 2-19 July 1990 to test 
the performance of the Naval Oceanic Vertical Aerosol Model (NOVAM) in a weak- 
convective regime. Meteorological data collected by aircraft and boat was used to gen- 
erate the surface and vertical profile information files required by NOVAM. Using this 
information, NOVAM predicts the aerosol extinction {knr^) for a vertical cross-section 
of the atmosphere. Aircraft-observed aerosol extinction profiles were also obtained. 
Comparisons between observed and NOVA.M aerosol extinction profiles revealed major 
deviations above the cloud top. From the surface to the top of the cloud layer. NOVA.M 
generally did an excellent job in predicting profile shape, with the magnitude of aerosol 
extinction tied to the extinction matched at the surface. In a few cases, observed ex- 
tinction increased more rapidly than NOVAM predicted extinction from the surface to 
the base of the cloud layer. This is attributed to rain scavenging associated with 
thunderstorm activity in the area. Comparison between different aerosol extinction 
profiles revealed much spatial and temporal variation that was verified by Lidar profiles 
of atmospheric structure. Thunderstorm activity, multiple cloud-layers, and the spatial 
variation in the atmospheric structure have led to a hypothesis that deep-convection was 
responsible for the major differences between observed and predicted aerosol extinction 
profiles. If this is the case, a simple modification to the weak-convective model may be 
made to apply to a deep-convection model. 



c m5 Ws- 

i./ 

TABLE OF CONTENTS 

I. INTRODUCTION I 

II. BACKGROUND 4 

A. EXTINCTION PARA .METERS 4 

B. THE LOW RESOLUTION TRANSMITTANCE MODEL (LOWTRAN) . . 5 

C. THE NAVY AEROSOL MODEL (NAM) 5 

D. THE MARITIME ATMOSPHERIC BOUNDARY LAYER (MABL) 7 

HI. THE NAVAL OCEANIC VERTICAL AEROSOL MODEL (NOVAM) 9 

A. BACKGROUND 9 

B. NOVAM INPUTS 10 

C. SUB-.MODELS 12 

I. The default profile generator 13 

IV. INFLUENCES ON AEROSOL SIZE DISTRIBUTION PROFILES WITH 

WEAK CONVECTION 15 

A. CONCENTR.ATION GR.ADIENTS IN NOVAM 15 

B. RELATIVE HUMIDITY EFFECTS IN NOVAM 18 

V. THE KEY90 EXPERI.MENT 21 

A. GOAL 21 

B. .MEASURE.MENTS 21 

C. ATMOSPHERIC SYNOPTIC SCALE 23 

VI. NOVAM USAGE, PROCEDURE AND RESULTS 27 

A. PROCEDURE 27 

1. Surface File Generation 27 

2. Radiosonde Preamble Generation 27 

3. NOVAM Runs 27 

B. RESULTS 28 



IV 



VII. SUMMARY, CONCLUSION AND RnCOMMCNDATlONS 47 

LIST OF REFERENCLS 49 

INITIAL DISTRIBUTION LIST 53 



LIST OF FIGURES 



Figure 1. KEY90 experimental area around Marathon (shaded area) 3 

Figure 2. Atmospheric Transmittance (0 - 15 microns) 5 

Figure 3. Idealized three lognormal marine aerosol 7 

Figure 4. Schematic of dominant processes in a multi-layered MABL 8 

Figure 5. Stylized profile definitions for use in preamble 13 

Figure 6. NOVAM flow diagram 14 

Figure 7. Schematic diagram of aerosol concentration 17 

Figure 8. Growth curves for marine aerosols 19 

Figure 9. Summary of instrumentation used during KEY90 22 

Figure 10. Synoptic pattern in Marathon area 24 

Figure 11. Synoptic pattern in Marathon area (cont.) 25 

Figure 12. Flight extinction profiles with flight composite 30 

Figure 13. NOVA.M prediction profiles with NOVAM composite 31 

Figure 14. 9 July results 33 

Figure 15. 10 July results 34 

Figure 16. 12 July results for 0030 UTC 35 

Figure 17. 12 July results for 1930 UTC (a) 36 

Figure 18. 12 July results for 1930 UTC (b) 37 

Figure 19. 13 July results 38 

Figure 20. 14 July results 39 

Figure 21. 16 July results for 2100 UTC 40 

Figure 22. 16 July results for 2359 UTC 41 

Figure 23. 17 July results for 1605 UTC 42 

Figure 24. 18 July results for 1240 UTC 43 

Figure 25. 18 July results for 1545 UTC 44 

Figure 26. 19 July results 45 

Figure 27. Logarithmically amplified XRL Lidar extinction profiles 46 



VI 



LIST OF TABLES 



Table 1. SURFACE OBSERVATION DATA FILE 10 

Table 2. CONTENTS OF THE RADIOSONDE DATA FILE PREAMBLE ... 1 1 

Table 3. GROWTH FACTOR CONSTANTS AND VALIDFI Y IGXNGE 20 

Table 4. SOURCES FOR VERTICAL PROl ILLS 23 



vii 



ACKNOWLEDGEMENTS 



As much as 1 would like to say 1 did it all myself, that simply is not the truth. At 
the top of my list, 1 must thank Professor Kenneth L. Davidson for providing his brand 
of eternal optimism, even when this thesis was in its darkest hours. Additionally, inputs 
from Professor (and fellow Seaslug) Philip A. Durkee and .Mathieu Rouault (when he 
wasn't out surfing!) proved to be invaluable. The most rewarding experience for me was 
meeting and working with the scientists involved in the KEY90 experiment. Special 
thanks to Stu Gathman and Gerrit de Leeuw for their personal contributions, in addition 
to the numerous references they provided. Obviously, without data provided by the 
KEY90 team, there wouldn't be many results to talk about, therefore I greatly appreci- 
ate the efforts of Doug Jensen, Bill Hooper, Keith Jones, and the many others involved. 
Last, but certainly not least, I wish to thank my family and friends that have supported 
me throughout. 



viii 



I. INTRODUCTION 



Successful employment of modern weapon systems require the continuous ability to 
use wavelengths in the electro-magnetic (EM) spectrum. This is the fundamental con- 
cept of electronic warfare and such importance has been placed on this mission area that 
the following ma.xim was developed; The side that controls the electro-magnetic (EM) 
spectrum controls the outcome of any conilict in modern war or global politics'. Oper- 
ation Desert Storm was successful for many reasons, one of which was the allied forces 
superiority in the maintained use of systems operating at various wavelengths through- 
out the EM spectrum. Even with this superiority, there is much room for improvement. 
One example is the limitation that cloud covxr had on operations. Another is operating 
in regions with extensive aerosol content (such as were caused by the intentional burning 
of oil rigs). 

Many military devices operate in the visible to the far infrared (IR) wavelengths. 
Some examples are forward-looking infrared (FLIR), precision guided munitions and 
other electro-optical (EO) systems. For this reason, there is a great deal of interest in 
reliable estimation of EO instrumentation performance for vertical and slant-path ob- 
servations. For systems such as those above, this requires knowledge of the vertical 
variation of aerosol scattering and absorption in the wavelengths of operation. 

LOWTRAN 7 (Low resolution transmittance model version 7) calculates atmo- 
spheric scattering and absorption caused by aerosols and molecules along a non- 
homogeneous path. Because there isn't one aerosol model that predicts scattering and 
absorption the best under all conditions, there are dilTerent aerosol models a user of 
LOWTR/\N may choose. For maritime applications, the present model is the Na\y 
Aerosol .Model (NAM). NAM attempts to relate aerosol size distributions to 
meteorological parameters such as relative humidity, wind speed, and visibility, and then 
calculate the optical properties of the modeled aerosol. Most empirical models, such as 
NAM, do not predict the vertical distribution of aerosol required to extend extinction 
predictions to higher levels. For the maritime regime, this led to the development of the 
Naval Oceanic Vertical Aerosol Model (N'OVA.M). 

The purpose of NOVAM (Gathman, 1989) is to determine the vertical distribution 
of the aerosol size spectra in the marine environment, used to obtain the associated op- 
tical and infrared properties along slant paths within the marine boundary layer. The 



means used to achieve this goal are to integrate various marine aerosol studies of a di- 
verse nature into a single model. This integration has led to the development of sub- 
models within NOVAM, one of which is selected based on the parameters that arc input. 

This thesis will investigate the performance of NOVAM under conditions of weak 
cumulus convection. These conditions existed during the KEY90 experiment which took 
place from 2 July 1990 through 19 July 1990, near Marathon, Florida on the Florida 
Keys (Figure 1). 
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Figure I. KEY90 experimental area around Marathon (shaded area) 
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II. BACKGROUND 



A. EXTINCTION PARAMETERS 

In this section, fundamental concepts and terms that form the basis of the 
LOWTRAN model, and led to the development of NAM and NOVAM are presented. 
Lambert, Beer and Bouguet established a relationship between attenuated radiation and 
the incident radiation for a homogeneous extinction as follow's: 

1= I^e\p{-b.x) (1) 



where 1 is the attenuated radiation, 

/„ is the incident radiation, 
b is the total extinction coelTicient, 

X is the path length. 

To determine the total extinction coefficient, b, knowledge that extinction is caused by 
absorption and scattering by molecules and aerosol (dry or liquid) contributing linearly 
is applied as follows: 

b — bf,j^ + b^^ + bg^ + bg^ ( 2 ) 

( 

where b„, = extinction due to molecular absorption, 

b,„, = extinction due to molecular (Rayleigh) scattering, 
b^^ = extinction due to aerosol absorption. 

= extinction due to aerosol (Mie) scattering. 

In the visible region of the spectrum, b^, and b„, are the dominant extinction com- 
ponents. In the IR regions, b„^ and are the primar>’ components causing lowered 
transmittance (the ratio of I to /„). Figure 2 on page 5 shows the effect molecular ab- 
sorption has on transmittance over a horizontal 1.82 km path for wavelengths up to 15 
nm. This figure clearly demonstrates that atmospheric total e.xtinction is w'avelength 
dependant due to at least its b„^ component. 

Equation 1 assumes a homogeneous total extinction throughout the atmosphere. 
Atmospheric studies have show'n this assumption to generally be valid for horizontal 
paths (such as Figure 2), however, a model needed to be developed to predict the vari- 
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Figure 2. Atmospheric Transmittance (U - 15 microns): Atmospheric 

Transmittance measured over a 1.82 km horizontal path at sea level 
(Hudson, 1969). 

ations of e.xtinction in the vertical required to perform slant-path calculations. 
LOWTRAN is the model that resulted from this requirement. LOWTRAN was de- 
signed to calculate the total e.xtinction as a function of altitude, and produce as output 
the radiance and/or the transmittance for a specified path through the atmosphere. 

B. THE LOW RESOLUTION TRANSMITTANCE MODEL (LOWTRAN) 
LOWTRAN's radiance calculations account for contributions from atmospheric 

self-emission, solar and/or lunar radiance single scattered into the path, direct solar 
irradiance through a slant path to space and multiple scattered solar and/or self- 
emmission radiance into a path. LOW'l'RAN 7 (the latest version) incorporates separate 
molecular profiles for all major, as well as 13 minor and trace gases found in the at- 
mosphere. In addition, various aerosol, cloud and rain models are included to allow the 
user to select the models that most closely represent the region of interest, fhe user also 
has the option to select from six reference atmospheric profiles, or define a new atmo- 
spheric profile. A more extensive description of LOW TRAN may be found in ONTAR 
corp. (1990). (Kneizys et al., 1988) 

C. THE NAVY AEROSOL MODEL (NAM) 

NA.M represents one of the several empirically derived models contained in 
LOW 1 RAN and was designed to estimate the contribution of aerosol to the LO prop- 
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agating characteristics in the marine environment. This model predicts an aerosol size 
distribution (dN dr) within the atmospheric sin face-layer, nominally for 10 meters, for 
specified wind speed (current and 24-hour averaged), visibility, and relative humidity 
(Rll) values. Three distinct types (modes) of aerosols are assumed to exist within the 
size distribution: 

1. A continental component contributed by a nearby land mass, 

2. A stationary component affected by the winds (controlled by the 24-hour averaged 
wind speed), 

3. A fresh component caused by the current wind over the water. 

The predicted aerosol size distribution is obtained by summing the contribution of each 
component as follows: 



where X represents the number density per unit radius of particles at radius ,r , 

A, is the amplitude parameter for mode i, 

/ is the aerosol growth parameter for mode i (to be discussed in a later section), 
is the mode radius parameter. 

When a standard relative humidity of S0% is introduced into the above equation, / is 
replaced by the value one, as will be shown in the relative humidity discussion. A plot 
of log(dX dr) at 80% relative humidity versus log(r) is shown in Figure 3. In this figure, 
.'\1, A2, and A3 are the amplitudesparameters as functions of air mass type, 24-hour 
averaged wind speed, and current wind speed respectively; and r is the aerosol particle 
radius at the ambient relative humidity. The mode radius parameter, r,„, is the radius that 
the maximum amplitude parameter value for each mode, i, occurs for this 80% relative 
humidity. (Gathman, 1983) 

NA.VFs application is to make horizontal path calculations of aerosol extinction at 
the sea surface based on the following relationship between dX/dr and aerosol ex- 
tinction (/?„,): 




(3) 



f^exi ^aa "b Qext ' 




(4) 
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Figure 3. Idenlized tliree lognormal marine aerosol: (Gathman, 1983). 

where Q,,, is the Mie eniciency coelficient for extinction which is a function of wave- 
length, causing /?„, to be wavelength dependent. 

To use NAM in slant path calculations, an assumed vertical size distribution must 
be developed for each mode based on the size distribution NAM produces at 10 meters. 
The present assumed size distribution profile has an exponential decrease with height. 
NOVAM is being developed to provide a capability for more accurate slant-path aer- 
osol extinction calculations in the maritime regime. 

D. THE MARITIME ATMOSPHERIC BOUNDARY LAYER (MABL) 

Vertical profiles of aerosol depend on the thermodynamic structure and processes 
within the MABL. A schematic representation of the MABL is shown in Figure 4. In 
this scheme, the cloud layer varies the most and, at times, may not be present. In cases 
where the large scale flow prevents cloud development, the MABL terminates at the 
transition layer. In cases where very deep cloud development has taken place, the top 
of the MABL is not well defined. F'eatures such as frontogenesis and horizontal 
advection would cause the one-dimensional view of the MABL {as depicted in 
Figure 4) to be no longer valid. 

For a MABL as shown in Figure 4, the Ibllowing principles apply: (Augstein, 1976) 
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Figure 4. Schematic of dominant processes in a multi-layered MABL: 0 = poten- 

tial temperature. (Augsteiii, 1976) 



1. Radiative cooling is compensated for mainly by the adiabatic warming of large- 
scale subsiding air and, to a lesser extent, by the sensible heat input from the ocean. 

2. Buoyancy driven mixing produces a well-mixed region above the shallow surface 
layer. 

3. Kinetic energy which is generated by the concurrent heating and moistening of air 
from below is partly converted to potential energy through entrainment of poten- 
tially warmer air from above through the top of the mixed layer. This assumes that 
buoyant parcels overshoot their level of equilibrium, thereby forcing mass of less 
density downward. 

4. The top of the mixed layer is determined by a balance of small-scale mass 
entrainment from below and the large-scale downward mass flow. 

5. Clouds begin to develop when convective plumes or bubbles penetrate to their 
condensation level. 
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III. THE NAVAL OCEANIC VERTICAL AEROSOL MODEL (NOVAM) 
A. BACKGROUND 

As described previously, NAM does not predict the correct vertical distributions 
which are required to extend extinction predictions to higher levels. NOVA.M is designed 
to predict non-uniform and non-logarithmic extinction profiles observed in the MABL 
for wavelengths between .2 and 40 nm. (Gathman, 1989) 

NOVAM is a combination of empirical and physical models which describe aerosol 
dynamical behavior. The current version of NOVAM incorporates the following mod- 
els; 

1. N.AM, which defines the surface (10 meter) aerosol characteristics (Gathman, 
1983), 

2. Micro-meteorological models, which describe the vertical distribution of aerosol in 
the well-mixed .MABL (Fairall and Davidson, 1986), 

3. A model that describes the vertical distribution of aerosol in the weak cumulus 
convection regime (described below) (Davidson and Fairall, 1986), 

4. A model describing the aerosol structure below stratus decks (Noonkester, 1985), 

5. A default boundary layer structure model of meteorological parameters based on 
surface observations (Gathman, 1989), 

6. A model that predicts aerosol growth based on relative humidity and the air mass 
type in the well-mixed and cloud layers (1 itzgerald, 1978; Gerber, 1985). 

Extinction properties are calculated from the generated aerosol profiles using .MIE 
scattering for which two basic characteristics exist: 

1. The aerosol particle radius and the wavelength of incident energy are of the same 
order of magnitude. 

2. Scattering is generally concentrated in the forward direction (forward scatter). 

The following particle source and particle physics principles have been used in the 
development of NOVAM. Large sea-salt aerosol particles are produced from white 
water phenomena at the air-sea interface (.Monahan et al., 1982). Smaller aerosol par- 
ticles are produced from sources such as gas-to-particle conversion or anthropologically 
generated aerosol. The concentration of any particular aerosol is dependant on the 
source strength of aerosol production and on the mixing process as it relates to scalar 
contaminants. Regarding hygroscopic sea-salt aerosol, its size is very dependant on the 
relative humidity it is immersed in. Also, as the hygroscopic aerosol picks up water vapor 
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from the atmosphere and grows in size, it changes its chemical composition and its index 
of refraction. 

B. NOVAM INPUTS 

Tables 1 and 2 show surface and radiosonde observation files (inputs) required by 
NOVAM. The surface observed parameters are used to perform the size distribution and 
aerosol extinction calculations at the surface and for sub-model selection, in conjunction 
with some information contained in Table 2. The cloud type, present weather, and the 
zonal seasonal category codes may all be found in separate tables as described by 
Gathman (1989). If any of these data items are unavailable for input, the model calcu- 
lates that parameter using default information which results in a reduced confidence in 
aerosol extinction predictions. 



Table 1. SURFACE OBSERVATION DATA FILE 



Position 


.Vleteorological Data 


1. 


Sea surface temperature (°C) 




Air temperature ("O 


3. 


Relative humidity (RH) f/o) 


4. 


Optical visibility (km) 


5. 


Current real wind speed (m s) 


6. 


Averaged wind speed (24 h) (m s) 


7. 


Air mass parameter (see text) 


8. 


Cloud cover (tenths) 


9. 


Cloud type [0...9] 


10. 


Surface IR ext. (I 'km) (,'t^ IC.b^t/n) 


11. 


Present weather in standard code [0...99] 


12. 


Height of lowest cloud (m) 


13. 


Zonal/seasonal categorv’ [1...6] 



Most of the parameters listed in Table I are self explanatory, however, a couple 
warrant further examination. The air mass parameter is an indicator of the degree of 
continental contaminants and may be represented as follows: 

a.m.p. = RnjA A- 1 

where Rn represents the atmospheric radon content (picocuries per cubic cm). 



10 



Table 2. CONTENTS OF THE RADIOSONDE DATA FILE PREAMBLE 



Row Column Description 


[ 1 ]‘ 


[ 1 ] 


Number of actual radiosonde observations 


[ 1 ]‘ 


[ 2 ] 


Surface potential temperature from radiosonde 


[ 1 ]' 


[3] 


Surface mi.xing ratio from radiosonde 


Profile Characteristics 


[ 2 ] 


[ 1 ] 


Height of base of cloud layer (CB) 


[ 2 ] 


[ 2 ] 


Potential temperature just below CB 


[ 2 ] 


[3] 


.Vlixing ratio just below CB 


[3] 


[ 1 ] 


Units indicator for temperature ( = 1 if “C) 


[3] 


[ 2 ] 


Potential temperature just above CB 


[3] 


[3] 


Mi.xing ratio just above CB 


[4] 


[ 1 ] 


Height of the cloud layer top (CT) 


[4] 


[ 2 ] 


Potential temperature just below CT 


[4] 


[3] 


.Mixing ratio just below C’f 


[5] 


[ 1 ] 


Units indicator for ntixing ratio ( = 1 ifg kg) 


[5] 


[ 2 ] 


Potential Temperature just above CT 


[5] 


[3J 


Mi.xing ratio just above CT 


‘ .Minimum requirement for radiosonde data file. 



or as. 



a.m.p. = (9 X exp[ — f/4]) + 1 

where t equals the elapsed time (in days) it takes the current air mass to reach the point 
of observation from a distant land mass. 

An air mass parameter of 1 indicates a pure air mass (no contaniinants). When the air 
mass parameter is less than or equal to 5, the aerosol is assumed to be made up of three 
lognormal components (/I,, A 2 , and /I 3 ) as contained in the NAM description above. 
For air mass parameters greater than 5, indicating the presence of non-soluble aerosol, 
an additional lognormal component (/Iq) is assumed. This class of non-soluble aerosol 
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is only present close to shore lines and has the same mode radius as shown in 
Figure 3. Being assumed non-soluble, the mode 0 component of aerosol does not change 
size with changes in relative humidity. (Gathman. 1989) 

The remaining inputs to NOVAM are contained in a data file describing the vertical 
atmosphere, referred to as the radiosonde data file in Gathman, 1989. Table 2 shows the 
format for the first five rows that are contained in this file, necessary when choosing the 
weak convective model. Figure 5 graphically illustrates the position in the file by row 
and column designation. The remainder of the radiosonde file consists of radiosonde 
observation height, potential temperature, and mixing ratio entries. If radiosonde data 
is not available, a minimum of 1 surface entry must be input and a default profile is se- 
lected. 

After the input parameters are translated by NOVAM to predictions of mixing, 
source strengths, and size distributions, a prediction is provided of the extinction of EO 
energy as a function of altitude for the wavelength of interest. This output is based on 
the dry aerosol size distribution at a particular height and adjusted for relative humidity 
at that height. 

C. SUB-MODELS 

This section is intended to give a brief explanation of the different models being de- 
veloped within NOVAM. Figure 6 shows the major decision points and the track of 
flow of information in the current version. There are currently four sub-models in 
NOVAM: 

1. A mixed boundary layer model (Fairall & Davidson, 1986; Davidson & Fairall, 
1986) where the boundary layer depth must be less than 3 km and not stable, 

2. A stratus model for winds less than 5 m s model (Noonkester, 1985). Here, the size 
distribution is not represented by a series of lognormals (as in Figure 3), there is 
no requirement for calculations to fit NAM predictions at 10 meters, and the 
wavelengths available for NOVAM calculation are from 1 to 1 1 The limita- 
tion on wavelengths is due to optical calculations using a simplification of Mie 
scattering parameters. 

3. A weak convection model (Davidson & Fairall, 1986) used when scattered cumulus 
clouds are present, cloud tops do not exceed 3 km, and a well-mixed layer is present 
below the cloud base, 

4. A default profile (Gathman, 1989) which is used when there is no vertical potential 
temperature and mi.xing ratio information. 

A stratus model for winds greater than 5 m/s and a deep convection model for cloud, tops 
extending above 3000 m are not yet supported. 



12 




Figure 5. Stylized profile definitions for use in preainhle: (Gatliinan, 1989). 

NOVAM sub-model selection is based on the input parameters describing the ver- 
tical stratification (thermal stability, presence of an inversion, and the inversion height), 
cloud cover, cloud type, wind speed, and the requested wavelength for extinction calcu- 
lations. 

I. The default profile generator 

The default profile generator uses the surface meteorological measurements to 
develop profiles of temperature and relative humidity, generating values for every 100 
meters of altitude. The default profile differs from a radiosonde derived profile by having 
the vertical distribution of the lognormal components described by an exponential 
function of scale height. For each default altitude level, the optical calculations are based 
on the combined size distributions of the aerosol deduced from the model assumptions. 
The index of refraction and size arc determined by the response (swelling) of hygroscopic 
aerosols based on the relative humidity calculated from the relative humidity profile 
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Figure 6. NOVAM flow diagrnin: (Gathnian, 1989). 



generator (Gathman, 78). From these parameters, the aerosol extinction at each level 
can be calculated directly. 
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IV. INFLUENCES ON AEROSOL SIZE DISTRIBUTION PROFILES 

WITH WEAK CONVECTION 



A. CONCENTRATION GRADIENTS IN NOVAM 

The basis for modeling the evolution of aerosol concentration is the continuity 

equation. In the Boussinesq approximation, which assumes a shallow planetary 

dp . . - 

boundary layer, an incompressible fluid = 0), and negligible horizontal gradients of 

density, an aerosol component concentration, X, can be written in its simplest one- 
dimensional form as 




-f L V^^A -f- IV 



cA 




-^((I'A'- V.X-D' 



G2 



where U is the mean horizontal wind, 

W is the mean vertical wind, 

V, is the particle gravitational settling velocity, 

D is the particle molecular transport coeiricient, 

li'X is the mean vertical flux of particles due to turbulent transport. 



( 7 ) 



In this equation, X must be a conservative aerosol variable, e.g. dry size concentration, 
that is unalfected by variations in relative humidity. Terms such as iVX are unlikely to 
be known from first principles leading to applying models of the boundary layer to sim- 
plify the process. (Davidson and Fairall, 1986) 

The Albrecht model was developed specifically for a trade-wind region (Albrecht, 
1979) and has been incorporated in NOVAM to describe the weak convective regime. 
The three-layer structure of this regime, described by Albrecht, is shown in Figure 7. 
In Figure 7, C represents the dry aerosol concentration with the subscripts r, b, c, and 
i representing the surface, cloud base, the in-cloud layer, and the cloud top respectively. 
The layers and their characteristics are summarized below, Davidson & Fairall (1986): 

1. The subcloud layer (0 < z < z*) where z^ is the cloud base (LCL). This layer is 
characterized by linear flux profiles and "well-ntixed'' mean profiles. The depth of 
this layer is changed by subsidence and entrainment. 

2. The cumulus cloud layer (z* < z < z,). characterized by parabolic fiux profiles and 
linear mean profiles. The slopes of the mean properties are determined by empirical 
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models of cloud transport and entrainment properties as well as the boundary 
conditions at the top and bottom of the cloud layer. 

3. The free troposphere (non-turbulent) above the capping "trade" inversion. The 
properties of this layer are determined by advection and are usually input variables 
to the model. 



Note that the mean profiles are not continuous at interfaces z* and z,. 

Factors that influence the in-cloud dry aerosol concentration gradient, y (shown in 
Figure 7), are given by the following balance equation: 



Dy 

Dt 



= (V.F- f )y -b 



E^X, 



-b 



4w',AA; 



(z,-z,) 



( 8 ) 



where V V is divergence. 

iv', is the entrainment rate of the inversion, 

T is the cloud relaxation time in days, 

E is the cloud wall entrainment variable given by; 

^ (2,-2,);v+ 2A7;(1 +n 



where AT„ = .S'K, 

F^l.6 , 

y, = the entrainment rate of B, in the layer, 
is the jump in 0, at the cloud base. 



Because the layer below the cloud base is considered well-mixed, the dry aerosol con- 
centration gradient for this layer is normally considered to be 0, but a non-zero gradient 
could be predicted by the model with the specified parameters. (Davidson and Fairall, 
1986) 

Some assumptions were made to simplify calculations when incorporating the 
Albrecht model into N'OVAM. The Albrecht model is dynamic, predicting the temporal 
evolution of the boundary layer. N’OVAM seeks to characterize the average or typical 
vertical structure, therefore all time derivatives are neglected in the above equations. 
Additionally, the following assumptions are used to simplify calculations; (Davidson & 
Fairall, 1986) 

1. In the mixed layer, C = C, 
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Figure 7. Schematic diagram of aerosol concentration: Dry aerosol concentration 

based on maritime and continental sources for a weak cumulus con- 
vection boundary layer structure. (Gathman, 1989) 



2. In the cloud layer, C = Q + — zj 

3. C is given by a standard wind speed dependent model 

4. At the cloud top, assume: 

a. C* = 0 if C represents a locally generated sea-salt mode 

b. C; is given by climatology 

5. Q = C; = Q 

These assumptions lead to the following equation for the vertical gradient from 
cloud base to cloud top: 



(cr~c> 

(1 +/t)Az 



( 10 ) 



where n is a dimensionless entrainment parameter at the cloud top given by: 
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n 



4w',x (r-y.r) 
Az 



( 11 ) 



In this form, the entrainment rate at the cloud base and the cloud wall entrainment are 
not important. These simplifications may reduce the accuracy in the vertical represen- 
tation of dry aerosol concentration, but this predicted vertical structure is a large first 
step in improving the exponential model currently in use. With mixing and entrainment 
controlled by the potential temperature and relative humidity profiles, a more accurate 
representation of the transport properties in the sub-cloud and cloud layers is obtained. 

B. RELATIVE HUMIDITY EFFECTS IN NOVAM 

Since aerosol particles attenuate energy in an amount proportional to their cross 
sectional area, any model requiring aerosol extinction calculations must take into ac- 
count the effect of relative humidity on aerosol growth. Vertical changes in relative 
humidity become important when requiring slant-path aerosol extinction calculations for 
the weak convective region where relative huntidity is expected to increase with height 
up to the cloud layer, then decrease gradually to top of the cloud layer. Figure 8 shows 
the efiect relative humidity has on particle radius (represented by the ratio of the aerosol 
particle radius at the indicated relative humidity (r) versus the radius of the same aerosol 
particle at 0% relative humidity (O). Fitzgerald found that pure NaCl particles experi- 
enced a sudden increase in size at some critical value of relative humidity (between 70 
and 76% relative humidity as represented by the hatched area in Figure S) whereas na- 
turally occurring maritime aerosol particles show continuous smooth growth with 
changing relative humidity. This may be explained by the gradual deliquescence of sol- 
uble material which is present as a mixture of different salts. 

Based on the information shown in Figure 8, Fitzgerald developed the following 
relationship: 



where y represents the air-mass characteristic and S represents the saturation ratio (for 
the purposes of this discussion, S can be assumed to be equal to RH / 100). 

The relative huntidity dependence in both NAM and NOVA.VI appears in the form 
of f, in the equation for size distribution (Equation 3). As can be seen in Figure 8. when 
the relative humidity increases, deviation of the measured ratios from any theoretical 




(12) 
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Figure 8. Growtli curves for marine aerosols: Theoretical and calculated growth 

curves for NaCl particles and natural marine aerosols (Fitzgerald, 1978). 



curve which is based on dry particle size becomes greater. Because relative humidity is 
gcncially above 50% throughout the MABL, a reference relative humidity of 80% was 
selected to minimize the error inherent in the theoretical aerosol growth calculation, and 
this reference relative humidity was used as a basis in computing/. An equation for/ 
is generated by dividing equation 12 by the ratio of r(.8) to r, yielding (Gerber, 1985); 



f, = 



r(S) 

/(. 8 ) 



C7,. - 5 1i 
C8,.(l - 5) ^ 



(13) 



where C7, C8 are conditional constants based on the aerosol mode component (recall 
Figure 3 on page 7) as displayed in Table 3. To incorporate the fourth mode (mode 0) 
into Equation 3, simply change the summation to include modes 0 through 3. 
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Table 3. GROWTH FACTOR CONSTANTS AND VALIDITY RANGE 



Aerosol Model 
Composition 


M ode 

Component 

(i) 


C7, 


C8, 


Range of 
Validity 


Sea salt 


3 


1.97 


5.83 


RH < 99.99% 


Sea salt 


2 


1.83 


5.13 


RH < 99 . 9 % 


Water Soluble 


1 


1.17 


1.87 


RH < 99% 


Dust 


0 


/o=1.0 


all values 
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V. THE KEY90 EXPERIMENT 



A. GOAL 

Most observational studies within the development of NOVA. M have been accom- 
plished oil' the California coast. Initial experiments involving the testing of NOVA.M 
were performed in the same environment that the model was developed in. As expected, 
the model verifies reasonably well in that environment. Since NOVAM is intended to 
be a global maritime model, the model had to be tested in other geographical scenarios. 

The KEY90 experiment was the first attempt to test NOVA .Vi's performance in a 
tropic-like environment where shallow or deep convection exists. The shaded region in 
Figure 1, known as .Marathon, was the region of data collection during KEY90. 
Insolation heating of the relatively shallow waters between Florida and Cuba causes 
convection to be the major meteorological process taking place near the air-sea interface 
in this region. Also, at .Vlarathon data could be obtained by both boat and aircraft while 
being away from land influences and major continental effects on the data. 

B. MEASUREMENTS 

Simultaneous measurements were made of parameters that the model requires as 
inputs and of extinction profiles at various wavelengths. The extinction profile meas- 
urements are then compared to the results obtained by running NOVA.M at wavelengths 
corresponding to those measured. The extinction profile measurements provided the 
standard to judge NOVAM extinction profile predictions. Although this procedure 
seems fairly straight forward, it is noted that it is difficult to provide an accurate estimate 
of either the measured extinction profile or of the quality of the meteorological data used 
as inputs to NOVA.M (e.g. surface IR extinction @ 10.6 nm and a.m.p.). 

The data measurements collected during KEY90 were made from shore, from a boat, 
and from two aircraft. Both the atmospheric surface layer and the .MADE structures 
were probed in detail. Figure 9 shows the instrumentation involved in the experiment. 
For each instrument, an asterisk marks the days of operation. Instruments are grouped 
according to the information provided (i.e. air temperature, sea surface temperature, 
relative hunridity, wind speed, wind direction, atmospheric pressure, visibility, vertical 
profiles, Lidar information, and other data). Table 4 gives a summary of the time and 
location that boat and aircraft atmospheric vertical profiles were taken. For easy refer- 
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Figure 9. Simiinary of instrunieiitation used during KEY90: (Gathman, 1991) 

ence, the difTerence of time (in minutes) and in location (in nautical miles) is provided 
in the last two columns. 
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Table 4. SOURCES FOR VERTICAL PROFILES 





.Aircraft 


Radiosonde 




Day 


Time 

(UTC) 


Location 


Time 

(U'fC) 


Location 


A t 
(min) 


A X 
(n.m. 


09 Jul 


2120 


A 


2200 


24°13'.V80°33TF 


40 


23 


10 Jul 


2100 


24'’39'A'S0°55' IF 


2200 


24"23'A' 80°33' IF 


60 


25 


12 Jul 


0030 


B 


2330' 


24'’2ILV 80"33'1F 


60 


26 


12 Jul 


1930 


B 


1906 


24M5'.V80°2r IF 


24 


39 




1945 


24°21'A' 80° 3 3' IF 


15 


39 


13 Jul 


1330 


B 


1319 


24°21 'A' 80°33' IF 


11 


26 


14 Jul 


1140 


B 


1 105 


24°2I'A' 80°33TF 


35 


26 


16 Jul 


2100 


A 


2U27 


24°13'.V 80°31'IF 




25 




2359 


B 


2324 


24°21'A'80°32' IF 


35 


27 


17 Jul 


1605 


B 


1617 


24°36LV 80°53' IF 


12 


2 


18 Jul 


1240 


B 


1218 


24°30'A' 80°53' IF 


22 


5 




1545 


A 


1459 


24°30'.V 8 l°O0' IF 


46 


9 


19 Jul 


1020 


B 


1023 


24°30LV80°53' IF 


3 


5 


’ 11 Jul 

A = 24“22'A'80°55'IF 
fl = 24“35';VS0°55'1F 



It should be noted that in the process of collecting data, areas of thunderstorm ac- 
tivity were avoided. This may have caused parameters such as average wind to be lower 
than what may have actually occurred. 

C. ATMOSPHERIC SYNOPTIC SCALE 

The major synoptic How patterns which influenced conditions in the .Marathon re- 
gion from 9 to 19 July are shown in Figure 10. For the period of 9 - 10 July, surface 
winds at Marathon were generally from the east to south-east with speeds of approxi- 
mately 5 m/s. This was a result of the influence of high pressure centered off the east 
coast of the U.S. (Figure lOa). The high pressure system continued to control the 
•Marathon region winds on the Ilth and I2th of July, however, short wave troughs were 
analyzed as passing through .Vlarathon on 11 July. 

As the high pressure system moved to the north-east, winds at .Marathon shifted 
from east-southeast to a more southerly direction (Figure 10b). Wind speeds were be- 
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Figure 10. Synoptic pattern in Marathon area: Surface pressure (mb) and 700 

mb heights (m * 10*‘ ) from NMC analysis for 1200 UTC is provided 
for tlie following days: a)9 July b)12 July. 
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Figure 11. Synoptic pattern in Maratlion area (cont.): Surface pressure (mb) and 

700 mb heights (m * 10 ‘ ) from NMC analysis for 1200 UTC is pro- 
vided for the following days: a)14 July b)17 July. 
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tween 3 and 5 m/s during this period. By 13 July, a low pressure system which was ap- 
proaching the eastern U.vS. had an influence on winds at Marathon. Wind flow remained 
southerly, but wind speed decreased slightly due to weaker pressure gradients over the 
area. On 14 July, the low pressure center had deepened and moved to the cast 
(Figure 11a). Surface wind flow over Marathon remained southerly with wind speeds 
increasing to 6 - 8 m/s. The cold front associated with the low pressure system was 
moving very slowly to the east. By late UTC on this day, it was evident that the frontal 
system was losing its upper level support, and the high pressure evident at 700 mb near 
the Florida Keys in the figure continued to build. 

From 15 to 19 July, high pressure once again dominated the flow pattern in the 
Marathon region (Figure 1 lb). As the high pressure center moved to the north-cast, the 
pressure gradient at all levels weakened. Surface winds were from the east to south-east 
throughout this period with wind speed decreasing to approximately 5 rn s by the 19th. 
A scries of short-wave troughs were analyzed as passing through the Marathon region 
on 17 July. 

Although all observations during KEY90 fell within the realm of weak convection 
defined by NOVAM, there were examples of deep convection in the form of 
thunderstorms. Thunderstorm activity was observed in the Marathon region on 14, 16, 
17, and 19 July. 
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VI. NOVAM USAGE, PROCEDURE AND RESULTS 



A. PROCEDURE 

1. Surface File Generation 

Each radiosonde and flight spiral profile requires corresponding surface file. 
Figure 9 shows that there were many choices for most surface file parameters. Because 
of consistent performance during the entire experiment, WeatherPak measurements were 
used for surface layer air temperatures, wind speeds, and relative humidities. When this 
information was not available on 10 and 1 1 July, sling psychrometer measurements were 
used for air temperature and relative humidity while a hand held device was used for 
wind speed. Bucket thermometer measurements of sea surface temperature were used. 
NRL's radon counter was used to establish the air mass parameter. The NOSC aircraft 
aerosol spectra was used to obtain the 10.6 fini aerosol extinction at the surface. The 
altitude of this reading was usually 30 meters above the surface. 1 he optical visibility 
value was calculated by the model, based on the above parameters. 

2. Radiosonde Preamble Generation 

The general shapes of the potential temperature and mixing ratio (for flight) or 
specific humidity (for radiosonde) profiles were exantined together when generating the 
radiosonde preamble. The first step was to establish heights of the first and second inv- 
ersions (zj and z, in Figure 5). In general, a well-mixed layer was easy to establish below 
the first inversion, where potential temperature and ntixing ratio remained constant. The 
more difiicult task was determining a height for the second inversion. In most cases, 
there was no clear boundary (jump) to define the top of the cloud layer, especially in the 
potential temperature profiles. As a result, relative humidity profiles and observations 
of slope changes in the potential temperature and mixing ratio profiles were considered 
subjectively in making the determination. 

3. NOVAM Runs 

NOVA VI runs for both the aircraft and radiosonde vertical profiles used the 
same surface file, so only one surface file was created for each row presented in 
Table 4. Information in the surface file is displayed above the radiosonde vertical pro- 
file panels in Figure 14 through Figure 26. 

The jumps occurring between straight line segments of the radiosonde and air- 
craft profiles are important run parameters. The predicted amount of entrainment oc- 
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curring at the cloud top (second inversion) is dependent on the difTerence (jumps) in the 
potential temperatures and mixing ratios at the level. The smaller the difference, the 
more entrainment is occurring between the cloud layer and the layer above the second 
inversion. The wavelength of interest for each model run was an input parameter since 
aerosol extinction varies with wavelength. All the results presented in this thesis are for 
a wavelength of 3.5 nm, even though aerosol extinction profiles were obtained by NOSC 
aircraft at .53 ixrn, 1.06 /xm, 3.5 fxm, and 10.6 txm. 

A NOVAM output file was generated for each vertical profile presented in Table 
4. A separate run was also made to generate a limiting profile. For each surface file, the 
NOVA.M default vertical profile was used with the relative humidity at each level set at 
98%. When aerosol extinction measurements in the lower portion of the MABL exceed 
this value, precipitation is believed to be a possible cause, and the assumption of 
NOVA.M that no precipitation exists is violated. The assumption of the default profile 
generator that each mode of the size distribution decreases exponentially with height 
will cause a difference in NOVAM versus limiting profile size distributions, increasing 
as the height of the well-mixed layer increases. 

The importance of the relative humidity to aerosol extinction can be shown by 
comparing NOVAM profiles to limiting profiles. At a wavelength of 3.5 ixm, an increase 
in relative humidity from the 80% typically observed at the surface to a 98% relative 
humidity (as used in the limiting profiles) results in an approximate order of magnitude 
increase in aerosol extinction prediction (comparison of aerosol extinction at surface 
between NOVA.M predictions based on vertical profiles and NOVA.M limiting profiles 
in Figure 14 through Figure 26). At the times of the vertical profiles, relative humidities 
above 90% were rarely seen, therefore a comparison could not be made of how well 
NOVA.M actually performed at the higher relative humidities. 

B. RESULTS 

Simple comparisons between the flight and radiosonde potential temperature, mix- 
ing ratio, and aerosol extinction profiles reveal extensive variability in space and time. 
Results from 19 July (Figure 26) illustrate this point the best. From Table 4, the ap- 
proximate horizontal distance between vertical profiles was five nautical miles, while the 
time difference was negligible. Analyses to determine the cloud-layer boundaries (defined 
by the first and second inversion) yielded a 150 meter difference in height of the cloud 
base and a 700 meter difference in the height that the second inversion occurred. The 
fact that there was a substantial spatial variability is further corroborated by the Lidar 
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extinction profile for 14 July shown in Figure 27, which clearly shows the spatial vari- 
ation in atmospheric structure. As discussed above, slant-path calculations of 
transmittance currently assume aerosol extinction to be horizontally uniform. This 
condition does not appear to exist in the KEY90 region. A predicted aerosol extinction 
profile generated by XOVAM based on one vertical profile will most likely not be rep- 
resentative of the average vertical structure over the path of interest. In order for 
XOVAM to provide the required predictions, either the mean cloud base and cloud top 
levels over the path need to be established, or a technique needs to be developed to 
minimize the importance of the analyzed inversion levels. 

All NOSC aircraft aerosol extinction observations between 9 - 19 July were plotted 
simultaneously (Figure 12) to examine how well XOVAM predictions characterized 
aerosol extinction. A subjectively determined mean profile was generated from the plots 
in Figure 12, and will be referred to as the flight composite. To be considered successful, 
separate XOVAM predictions should compare favorably to the llight composite under 
most circumstances. Due to space and time dilTerences between the flight profiles and 
the radiosonde profiles, only XOVAM results from flight profiles were used in the com- 
parisons with the observed aerosol extinction profiles. When considering only the top 
row of Figure 14 through Figure 26, the following features can be seen: 

1. XOVAM aerosol extinction predictions above the second inversion were extremely 
poor. 

2. Below the second inversion, XOVAM predictions compared favorably with the 
flight composite except in the cases where aerosol extinction dropped off noticeably 
(an order of magnitude difference between the surface and cloud base observations) 
below the cloud base (see top panel Figure 17c, Figure 20c, Figure 24c, and Fig- 
ure 25c). 

The assumption that aerosol extinction above the second inversion is essentially 
zero when the air mass parameter is calculated to be approximately equal to one does 
not appear to be a good one for this region. Aerosol extinction was determined to be 
significant if it is within one order of magnitude of the aerosol extinction that XOVAM 
calculates at the surface for the wavelength of interest. Using this as a criteria, the 12, 
13, 14, 17, and 18 July XOVA.M llight profiles all showed significant aerosol extinction 
contributions above the second inversion. 

One possible explanation for the poor performance of XOVA.M above the second 
inversion is that the weak or shallow convection assumptions are not elTective here ,i.e. 
this may be a case of deep convection. Thunderstorm activity and rain were frequently 
observed in .Marathon, and multiple cloud layers were not uncommon. The hypothesis 




is that relatively small regions of deep convective activity result in an increased pro- 
duction of aerosol at the surface and transports aerosol to higher levels. As these deep 
convective cells dissipate, large amounts of aerosol are still being transported horizon- 
tally throughout the region, which would account for the aerosol present above the 
second inversion in many cases. 

Another possible c.xplanation may be that the atmosphere above the second inver- 
sion commonly contains a significant number of small radii aerosol particles. These 
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Figure 13. NOVAM prediction proliles with iS'OVAM composite: 



particles may originate either from continental sources or from the ocean surface in the 
form of DMS. These DMS particles would be carried to the upper atmosphere by 
convective cells and transported by large scale circulation. Particle processes associated 
with DMS origins arc not well understood. (Durkee ct al., 1990) 

If the predicted aerosol e.xtinctions above the second inversion are ignored, and a 
simultaneous plot of NOVAM predictions based on flight vertical profiles is made, a 
NOVAM composite can be obtained (Idgurc 13) in the same manner as the Ilight 
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composite described above. The flight and NOVA.Vl composites turned out to be nearly 
identical (Figure 14c). There was not much variation in wind speed at Marathon and 
no frontal systems passed through during the KKY90 e.'cperiment, therefore it is likely 
that aerosol generation and transport mechanisms remaiined fairly constant during 
KEY90. This would e.xplain why composite mean profiles performed as well as NOVAM 
predictions below the second inversion. It is possible that mean profiles, such as the 
composites generated above, can characterize the typical aerosol extinction (at a 3.5 
nm wavelength) for the .Marathon region during July. This approach leads to the region 
being described empirically, therefore the predictions are specific to one region during 
one season. The required data base for describing all ocean regions for each season 
would be enormous and the time and expense involved would be substantial. 

'fhe fact that a simple exponential curve developed from composite mean profiles 
performed much better than NOVA.M above the evaluated second inversion indicates a 
need for further research on aerosol size distributions in the upper atmosphere. If the 
structure observed in KEY90 is in fact due to deep convective processes, a revision 
would have to be made to the sub-model selection process in NOVAM. Also, instead 
of the weak-convective assumption that the mode 2 and 3 components drop off to 0 
above the second inversion, an exponential function could be used to describe the aer- 
osol size distributions for these modes. 

It is believed that the observations where aerosol extinction decreased 'noticeably 
below the first inversion are attributed to the combination of increased aerosol pro- 
duction at a time prior to the observation, possibly caused by thunderstorm activity, and 
scavenging in the lower layer due to rain, creating a deficit of aerosol below the cloud 
layer. If future experiments reveal that rain scavenging accounts for reduced aerosol ex- 
tinction below the cloud layer, the present weather surface input file parameter is cur- 
rently the only means that this information can be relayed to NOVA.M. Because the 
aerosol extinction predicted at the surface is tied to surface observations of aerosol ex- 
tinction, the task of modeling this phenomenon would still remain. 



32 



UlXING RATIO (o/^i?) 




Figure 14. 9 July results: 'I he top panels are based on flight vertical profiles 

observed 09 July, 2120 U I'C. Bottom panels are based on radiosonde 
profiles for 09 July, 2200 U I C. Panels labeled a) represent temperature 
and mixing ratio (flight) or specific humidity (sonde) versus altitude, 
b) plots relative humidity versus altitude, and c) shows aerosol ex- 
tinction at 3.5 (im versus altitude. Refer to text for a description of the 
profiles appearing in panel c). 
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Figure 15. 10 July results: Tlie top panels are based on flight vertical profiles 

observed 10 July, 2100 UTC. Bottom panels are based on radiosonde 
profiles for 10 July, 2200 UTC. Panel description is same as in 
Figure 14. 
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Figiiie 16. 12 July results for 0030 UTC: Ihe top panels arc based on flight 

vertical profiles observed 12 July, 0030 LH C. Bottom panels arc based 
on radiosonde profiles for 11 July, 2330 U fC. Panel description is 
same as in Figure 14. 
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Figure 17. 12 July results for 1930 UTC (a): The top panels are based on flight 

vertical profiles observed 12 July, 1930 UTC. Bottom panels are based 
on radiosonde profiles for 12 Jul, 1906 UTC. Panel description is same 
as in Figure 14. 
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Figure 18. 12 July results for 1930 UTC (h): The top panels are based on flight 

vertical profiles observed 12 July, 1930 UTC. Bottom panels are based 
on radiosonde profiles for 12 .lul, 1945 UTC. Panel description is same 
as in Figure 14. 
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Figure 19. 13 July results: The top panels arc based on flight vertical profiles 

observed 13 July, 1330 UTC. nottoin panels arc based on radiosonde 
profiles for 13 July, 1319 UTC. Panel description is same as in 
Figure 14. 
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I 

Figure 20. 14 July results: The top panels are based on flight vertical profiles 

observed 14 July, 1140 U'l'C. Bottom panels are based on radiosonde 
profiles for 13 July, 1105 U'fC. Panel description is same as in 
Figure 14. 
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Figure 21. 16 July results for 2100 UTC: The top panels are based on flight 

vertical profiles observed 16 July, 2100 UTC. Bottom panels are based 
on radiosonde profiles for 16 July, 2027 UlC. Panel description is 
same as in Figure 14. 
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Figure 22. 16 July results for 2359 UTC: The top panels are based on flight 

vertical profiles observed 16 July, 2359 U'lC. Bottom panels are based 
on radiosonde profiles for 16 July, 2324 U’lC. Panel description is 
same as in Figure 14. 
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Figure 23. 17 July results for 1605 UTC: The top panels arc based on flight 

vertical profiles observed 17 July, 1605 UTC. Bottom panels arc based 
on radiosonde profiles for 17 July, 1617 UTC. Panel description is 
same as in Figure 14. 
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Figure 24. 18 July resuKs for 1240 UTC: The top panels arc based on flight 

vertical pronics observed 18 July, 1240 UTC. Bottom panels are based 
on radiosonde profiles for 18 July, 1218 UTC. ranel description is 
same as in Figure 14. 
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Figure 25. 18 July results for 1545 UTC: The top panels are based on flight 

vertical profiles observed 18 July, 1545 UTC. Bottom panels are based 
on radiosonde profiles for 18 July, 1459 UTC. Panel description is 
same as in Figure 14. 
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Figure 26. 19 July results: Tlic top panels arc based on (light vertical profiles 

observed 19 July, 1020 U’l C. Bottom panels arc based on radiosonde 
profiles for 19 July, 1023 UTC. Panel description is same as in 
Figure 14. 
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Figure 27. Logaritliinically ami)lified NRL Lidar extinction proliles: Coded in 

false color, horizontal axis is distance covered by NRL aircrall from 
1K)7 - 1126 U rC recorded on 14 July 1690. Vertical axis is altitude 
above sea level. Colors represent aerosol extinction in km ‘ as given by 
the bar code on the right hand side, (llooper, 1991) 
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VII. SUMMARY, CONCLUSION AND RECOMMENDATIONS 

An experiment was conducted in the Florida Keys from 2-19 July 1990 to test the 
performance of the Naval Oceanic Vertical Aerosol Model (NOVAM) in a tropical, 
weak-convective regime. Meteorological data collected by aircraft and boat was used to 
generate the surface and vertical profile information files required by NOVAM. Using 
this information, NOVAM predicts the aerosol extinction (Am~‘) for a vertical cross- 
section of the atmosphere. When NOVA.Vl is incorporated into the LOWTRAN model, 
which also accounts for molecular extinction, slant-path calculations for transmittance 
may be performed. 

Aircraft-observed aerosol extinction profiles were obtained in order to verify 
NOVA.M predictions of aerosol extinction. Comparisons between observed and 
NOV.AM aerosol extinction profiles revealed major deviations above the cloud top. 
From the surface to the top of the cloud layer, NOVA.M generally did an excellent job 
in predicting profile shape, with the magnitude of aerosol extinction tied to the e.x- 
tinction matched at the surface. In a few cases, observed extinction increased more rap- 
idly than NOVAM predicted extinction from the surface to the base of the cloud layer. 
This is attributed to rain scavenging associated with thunderstorm activity in the area. 
Comparison between different aerosol extinction profiles over time and space revealed 
much variation, which Lidar profiles of atmospheric structure verified. Due to these 
variations, many vertical soundings of the atmosphere would have to be obtained in or- 
der to characterize a mean atmosphere for the region of interest. Thunderstorm activity, 
multiple cloud-layers, and the spatial variation in the atmospheric structure have led to 
a hypothesis that deep-convection was responsible for the major differences between 
observed and predicted aerosol extinction profiles. If this is the case, a simple modifica- 
tion to the weak-convective model may be made to apply to a deep-convection model. 

The results described above, clearly show that NOVA.M is not yet ready to be used 
operationally in regions where convection is the dominant process. The following are 
recommendations for consideration in the improvement and development of NOVA.M: 

1. In order to eliminate the bias in wind speed measurements caused by avoiding areas 
of storm activity, wind measurements from buoys placed strategically in thd area 
of operations would be invaluable in generating the 24 hour average winds. 
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2. Conduct experiments in the MARA'l'HON area and other areas where weak 
convective conditions are expected. The data and results from these experiments 
could then be compared with the data and results from the KEY90 experiment. 

3. Continue to explore the possibility that some cases that are currently classified as 
weak convection may in fact be examples of deep convection. If this hypothesis is 
correct, the parameters used to difi'erentiate between shallow and deep convection 
need to be modified. 

4. Explore the possibility of using satellite information to improve model predictions. 
One example would be to use the techniques established by Durkee, et al. (1986; 
1990) to determine variations in optical depth (total extinction) over a horizontal 
path which could then be incorporated into NOVA.M. 

5. A comparison needs to be performed between what the current model is producing 
(LOWTFIAN N/\M) against what NOVA.M, incorporated in LOWTRAN, predicts 
(transmittance predictions versus measured would be extremely useful). 
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